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We prove that magic states from the Clifford hierarchy give optimal solutions for tasks involving 
nonlocality and entropic uncertainty with respect to Pauli measurements. For both the nonlocality 
and uncertainty tasks, stabilizer states are the worst possible pure states so our solutions have an 
operational interpretation as being highly non-stabilizer. The optimal strategy for a qudit version 
of the Clauser-Horne-Shimony-Holt (CHSH) game in prime dimensions is achieved by measuring 
maximally entangled states that are isomorphic to single-qudit magic states. These magic states 
have an appealingly simple form and our proof shows that they are “balanced” with respect to all but 
one of the mutually unbiased stabilizer bases. Of all equatorial qudit states, magic states minimize 
the average entropic uncertainties for collision entropy and also, for small prime dimensions, min- 
entropy - a fact that may have implications for cryptography. 


I. INTRODUCTION AND MOTIVATION 

In the context of fault-tolerant quantum computing us¬ 
ing error-correcting codes, there exist single-particle pure 
states -“magic states” - with desirable properties. Hav¬ 
ing access to a supply of these states expands the sets of 
robustly implementable operations to encompass a uni¬ 
versal gate set. Moreover, impure versions of these states 
can be purified using only the fault-tolerant operations 
provided by the error-correcting code. There is a privi¬ 
leged family of magic states, well-defined for both qubits 
and qudits of odd prime dimension p, whose structure 
reflects a group-theoretical object called the Clifford hi¬ 
erarchy. Here we will show that the utility of these magic 
states extends beyond the arena of quantum computa¬ 
tion; they essentially provide optimal strategies for both 
(i) a nonlocal game wherein both Alice and Bob each use 
p Pauli measurements, and (ii) a cryptographically mo¬ 
tivated scenario wherein we try to minimize the average 
entropic uncertainty with respect to a set of p mutually 
unbiased measurements. Qudits can have operational ad¬ 
vantages for nonlocal or cryptographic tasks in terms of 
tolerating inefficiencies or the addition of white noise [1- 
3], but the geometry of state space becomes notoriously 
complex when we venture past the qubit Bloch sphere. 
In this work we prove that the optimal states have a 
simple structure and furthermore that these states have 
connections to group-theoretical and number-theoretical 
structures. These ideas may prove useful elsewhere. 

Given the fundamental importance of the Clauser- 
Horne-Shimony-Holt [4] (CHSH) game in quantum infor¬ 
mation theory [5-7], it is well-motivated to study gener¬ 
alizations to different numbers of measurement settings, 
parties and so on. By preserving the structure of the 
CHSH game while enlarging the size of the input/output 
alphabet to d symbols (using qudits), we can examine 
how the quantum advantage scales with dimension or 
investigate new features that arise due to additional de¬ 
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grees of freedom. Buhrman and Massar [8] studied such a 
generalization and found bounds on the allowable quan¬ 
tum value. Ji et al. [9] constructed a Bell operator 
from qudit Pauli measurements, whose maximization de¬ 
scribes a quantum strategy for the qudit CHSH game. 
For small prime dimensions they analyzed classical (i.e. 
local hidden variable) and quantum values. A follow¬ 
up work by Liang et al. [10] expanded on this analysis, 
and broadened the numerical search to allow non-Pauli 
measurements. Most recently. Bavarian and Shor [11] 
looked at a generalized CHSH game wherein the alpha¬ 
bet consists of all field elements where g is a prime 
power, and proved a generalization of Tsirelon’s bound 
using Information Causality [12]. As was done by Ji et 
al. [9], we will restrict our analysis to projective Pauli 
measurements. Quantifying nonlocality is a tricky task, 
particularly if one is motivated by operational consider¬ 
ations. For example, an experimental test may achieve 
higher statistical significance via lower statistical error 
[13] or higher statistical strength via larger statistical di¬ 
vergence with any classical model [14]. One conceptually 
simple metric for nonlocality is to measure the amount 
by which a Bell inequality is violated. If Bell inequality 
violation is the only figure of merit then restricting to 
Pauli measurements is suboptimal; larger violations are 
possible with non-Pauli measurements whenever p > 3. 
Nevertheless, we remain particularly interested in non¬ 
stabilizer states, measurements (see e.g [15]) or transfor¬ 
mations that emerge as privileged with respect to their 
stabilizer counterparts, due to some geometrical or op¬ 
erational relationship. This is motivated in part by the 
clean mathematical and conceptual framework provided 
by the so-called “stabilizer subtheory” [16-18] of quantum 
mechanics for odd-dimensional qudits - loosely speak¬ 
ing, states and operations far outside this subtheory are 
highly non-classical [19]. 

Mutually unbiased bases (MUBs) are a crucial com¬ 
ponent in a variety of quantum-informational settings. 
They are necessary for the creation of maximally strong 
entropic uncertainty relations [20], which describe limits 
on the allowed probability distributions associated with 
multiple measurements on quantum states. In crypto- 
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graphic scenarios, uncertainty relations can be used to 
bound an eavesdropper’s knowledge of quantum states 
transmitted between two parties. In [21] Amburg et al. 
motivate the study of “MUB-balanced states” i.e., pure 
states for which, relative to a complete set of d+l MUBs, 
the list of probabilities of the d outcomes of one of these 
measurements is independent of the choice of measure¬ 
ment, up to permutations. Recent work by Appleby et 
al. [22] suggests that MUB-balanced states only exist in 
odd prime-power dimensions of the form d = 3 mod 4. 
A MUB-balanced state is automatically also a collision- 
entropic minimum uncertainty [22, 23] (or maximally cer¬ 
tain [24]) state, important in cryptography, but balanced¬ 
ness is a stricter condition than minimum uncertainty. 
From a foundational point of view, MUB-balanced states 
are analagous to harmonic oscillator eigenstates insofar 
as the “direction” of a measurement is unimportant. Our 
results include a concise proof that magic states are bal¬ 
anced with respect to p out of a total of p-l -1 total bases in 
a complete set of MUBs. This suggests that magic states 
may be relevant in quantum cryptography, in analogy 
with the Breidbart basis for qubits. We also comment on 
the connection between magic states and the Sato-Tate 
distribution that arises in number theory. 


II. BACKGROUND MATERIAL 
A. Prerequisites and Notation 


mutually unbiased bases (ii) the complete set of stabilizer 
states (Pauli eigenstates), for a single qudit. 

Consider the family of magic states / and magic gates 
M (diagonal in the computational basis) [25-27] 


I fa.b^c} 


= 4- p>3 


— O. G Z*, b,CG Zp 


( 4 ) 

( 5 ) 


where j -I-) is the equal-weighted superposition of all com¬ 
putational basis states and Z* = {1,2,... ,p — 1}. Note 
that the above expressions require a minor modification 
for the smallest prime dimensions [27] i.e.. 


\fa,b,c} — 



(| 0 ) -P7“+2f'+4c|l)) 
(|0) -p^2a-H6h-h3c|i^ 


a+66+6c 


| 2 )) 


p = 2 

p = 3 

( 6 ) 


where 7 = and ^ and a S Z*, 6 , c G Zp. 

For a fixed value a G Z* the magic states, in addition to 
the computational basis, form a complete set of MUBs 
since \{fa.b.c\fa.b'.c'}\ ~^b.h'^c,c' 4” (1 ^b.h')/\dp [28, 29]. 

Finally, define the Jamiolkowski state \Ja,b,c) as that 
which is created by applying a magic gate to one half 
of a Bell pair: 


The Weyl-Heisenberg operators T'(xlz) are a general¬ 
ized p-dimensional version of qubit Pauli operators 

Dix\z) = + (1) 

k 

with X\j) = \j + 1) Z\j) = w-^lj) (w = , 

where we always treat expressions like 4 to mean the 
mutiplicative inverse a~^ G Zp = {0,1,... ,p — 1}. The 
choice of phase has the convenient feature that D? , , = 

[X\Z) 

D{nx\nz) and we use symplectic notation to keep track of 
X and Z powers separately, so that multi-qudit operators 
look like 


D{xi\zi) Zl D(x2\z2) — D(^xi,X2\zi,Z2)- 

The rank -1 projector onto the eigenspace of D{i\b) 
(in other words, the projector onto the U-th vector of the 
B-th Weyl-Heisenberg basis) is given by 

= ( 2 ) 

3 

so that \iIjI) = — V (3) 

The p^ states V G Zp}, along with those of the 

computational basis, comprise both (i) a complete set of 


lJa,b,o} =^p<S>Ma,b,cl^}, 
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\jj) 

Vp' 


( 7 ) 


Clearly \Ja,b,c) is isomorphic to a single-qudit magic state 
|/a.h,c) under the identification \jj) G O \j) G C^'. 

The group of operations generated by (tensor products 
of) the Weyl-Heisenberg operators in Eq. (1) is the Pauli 
group (we denote it Ci for reasons that will become clear). 
The set of unitaries that map the Pauli group into itself 
under conjugation is the Clifford group (C 2 ). The uni¬ 
taries that map the Pauli group into the Clifford group, 
under conjugation, define [30] the third level of the Clif¬ 
ford hierarchy, C 3 , and so on. Overall we have 

Cfe := {U\UC,U^ C Cfe.!} ( 8 ) 


where Ci C C 2 C C 3 ... but the sets Cfc >3 no longer form 
a group. The operations in the first two levels arise nat¬ 
urally in the context of fault-tolerant quantum compu¬ 
tation, whereas an operation from third level is typically 
also required to enable universal quantum computation. 
Properties of C 3 for multiple qubits are discussed in [30] 
whereas the diagonal single-qudit subset of C 3 , as dis¬ 
cussed in [27], is given by Ma,b,c above. It is in this sense 
that we say states \fa,b,c) and | Ja,fc,c) reflect the structure 
of (the third level of) the Clifford hierarchy. 
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B. Generalizing the CHSH game 

The familiar two-qubit CHSH inequality takes the form 
('®)max — 2 with a Bell operator 

1 

B = XX + XY + YX -YY = ^ uj^^A^By (9) 

x,y—0 

where {Hq, Hi} = {Hq, Hi} = {X, F} ( 10 ) 

and (jj = — 1 is a p-th root of unity. The maximum ex¬ 
pectation value using local hidden variables (LHV) is 2 
whereas quantum mechanics can achieve 2-\/2 which is 
the maximum eigenvalue, A„iax(^), of the Bell operator. 
The optimal shared state that Alice and Bob can measure 
in this case is given by |Ti_o.o) from Eq. (7) i.e., 

(^) 2 “x=(^i.o.o|S|^i.o.o)= 2 y 2 . ( 11 ) 

Generalizing the above CHSH Bell operator to qudits 
[9], we hnd it convenient to define two closely related 
operators 

- E E , ( 12 ) 

^ nGZ* x,yGZp 

B=- (13) 

^ n,x,y^‘Ijp 

where w = so that B — B* Yplyi and consequently 

Amax(‘B) = Ai„ax(^*) + P- The traceless version B* is 
more convenient to work with but the maximal value of 
H is a more useful quantity. 

To aid our analysis we use the Pauli measurement op¬ 
erators adopted in [9] 

H, = By = (14) 

where x,y G Zp, so that Ay /2 = By. The operators 
G Zp} are a generalization of the qubit mea¬ 
surements {X,Y} used above insofar as they are the 
Pauli observables (excluding Z) whose eigenbases form 
mutually unbiased bases. In the qudit case, different or¬ 
derings and phases of the measurement operators can 
have an effect on the quantumly achievable expectation 
value. 

The expectation value of B from Eq. (13) can be related 
to the quantum value, 0 < ;/ < 1 , of a nonlocal game 
[5, 6 ], which takes the general form 

ly = ^ P{x,y)V{a,b,x,y)P{a,b\x,y) (15) 

a,b,x,yGZp 

for measurement settings {x,y), outcomes (a, 6 ) and a 
predicate or pay-off function V G {0,1}. By choos¬ 
ing a flat probability distribution over all measure¬ 
ment settings, P{x,y) = 1/p^, and a winning condition 
V (a, b, x, y) = 5a+h+xyp we arrive at 

i/= 4(^l^l^) = 4 E (16) 

a-\-h-\-xy—^ 


where the right hand side of the above expression is a 
more general form for CHSH-type games, amenable to 
POVMs as well as projective measurements (see also a 
discussion in [10]). 

The maximum classical value is given by find¬ 

ing local hidden variable assignments 

d — (uq , . . . , Qx 7 ■ ■ • j Oip— 1) € Zp 

b = (6 o,...,6p,...,6 p-i) G ZP 

that tell Alice and Bob which element of the spectrum of 
their observable they should output for a given measure¬ 
ment setting i.e., 

Ax I— y , By I— y . 

The best choices of a and b maximize the quantity 

n,x,y—0 

~ 'y ^a^+by+xy,0 

a,b.,x,y^'Ijp 

i.e., they give the best classical strategy for participants 
trying to maximize the number of instances where a -I- 
b + xy = 0 mod p. Liang et al. [10] provide an explicit 
classical strategy achieving 

(H)LHV>3p_2. 

Turning to quantum strategies. Bavarian and Shor’s 
result [11] (restricted to q = p) says that the quantum 
value (16) of a generalized CHSH game is bounded above 
(via Information Causality [12]) by 

" - '= K ' + gy) 

Clearly this implies {B)^^ < p'^IC{p) but it also ap¬ 
plies to more general scenarios wherein Alice and Bob 
use POVMs. It is also shown [11] that classical strate¬ 
gies can achieve H if q = or O ^ 

q = p^^~^ for a very small (< i^) positive constant 
Co- The moral is that we might not expect a classical- 
quantum separation in the asymptotic limit. Liang et 
al. [10] used the Navascues-Pironio-Acin [31] hierarchy 
to numerically bound the quantum value v from above 
and it appears that the values so obtained coincide with 
IC{p) for p = 5,7. 


C. Entropic Uncertainty Relations and 
Eavesdropping 

Fix the notation for the mutually unbiased basis 
(MUB) expansion of a generic operator K acting on a 
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qutrit state space as follows: 



C0,C5O 
Cl, 00 

Co,0 

Cl,0 

Co,l 

Cl,l 

Co,2 \ 
Cl.2 


y C2,cxD 

C2,0 

C2,l 

C2.2 / 


where 


cv,B = (I^b) from Eq. (3)) 

and the infinity basis label corresponds to the computa¬ 
tional basis i.e., IV’to) = 1^)- The generalization to higher 
dimensions is obvious. The coefficients cv,b correspond 
to probabilities when id is a normalized density operator, 
and the probability distribution associated with a basis 
B is given by {co,b, ci^s,..., Cp-i^s}- Inverting this de¬ 
composition, we find [32] 

B,V 

With the above decomposition in hand we can charac¬ 
terize a density matrix p in terms of its associated prob¬ 
ability distributions. The Renyi entropy of order 2 (i.e., 
the collision entropy) with respect to a particular basis B 
is given by H^ip) = — log {J2v^vb) whereas the min- 
entropy is given by = - log maxy cv,b [23, 24]. 

For an arbitrary pure state j^) G this can be rewritten 
as 



FIG. 1. The Breidbart basis {|/i,o,o): |/i,o,i)} (which is equal 
to {|R), jR^)}} in the magic state notation of Bravyi and 
Kitaev [52]) constitutes the optimal measurement basis for 
Eve if her goal is to infer as much information about Alice’s 
state as possible. The cryptographic protocol [37-39] assumes 
Alice is sending one of the four equatorial Pauli eigenstates. 


Now, instead of collision entropy, consider the total 
min-entropy across a number of unbiased measurement 
bases. The following min-entropic uncertainty relation 
was proven for dimensions of the form d = 2" in [24] and 
generally in [34] 


^ -id+l)\og 

B6{oo,Zd} 



d-1 Y 
Vd^)_ ■ 


(!</>)) = - log 

\ V 

^^fi„(l</>)) = -log maxl((/>lV'B)l^ 

which obey 

log(p) > > idfind^)) > 0. 

Typically we are interested in total or average entropies 
across a number of different measurement bases. The 
trade-offs that arise when we try to minimize entropy 
across multiple measurement bases gives insights into the 
structure of quantum state space. 

For a complete set of p-|-l unbiased measurement bases, 
the total collision entropy associated with any pure state 
is bounded from below by [23, 33] 

H -log5Il^‘^l^s)l'‘>-(T+l)log-^ (18) 

BGfoo.Zp} V ^ 

The factor of 2 in the numerator arises from the sur¬ 
prising fact that V 1"^ = 2 for any pure 

state {(j)). The above inequality is saturated whenever 
J2v Is independent of the basis B, as occurs 

for SIC-POVM fiducial states [33]. For this reason fidu¬ 
cial states qualify as “minimum-uncertainty states” (not 
all minimum uncertainty states are fiducial vectors, how¬ 
ever). MUB-balanced states [21] have the same probabil¬ 
ity distribution for all p -|- 1 bases and consequently also 
saturate this lower bound. 



Later we will interested in the case there we omit the 
computational basis and restrict to prime-dimensional 
systems so the relevant lower bound becomes 


Y. ^min(l<^)) > -plog 

BGl,p 



(19) 


As well as constraining quantum state space, the above 
entropic quantities have an operational meaning in the 
context of cryptography [20, 23, 35, 36]. As explained 
in e.g., [37], if the four possible signal states transmit¬ 
ted from Alice to Bob in the BB’84 quantum key distri¬ 
bution protocol [38] are the eigenstates of ax and ay, 
then the optimal measurement basis (see Fig. 1) for 
Eve’s intercept-resend attack is the basis of magic states 
{|/i,o,o)) l/i.o,i)} better known as {\H), \H^)} in the con¬ 
text of quantum computation or the Breidbart basis in 
cryptographic settings [37-39]. This basis, which can be 
thought of as splitting the difference between the ax and 
ay bases allows Eve to ascertain Alice’s information with 
probability 5(1 + ~ 85%, (which is the optimal value 

of the CHSH nonlocal game) and the state min¬ 

imizes the average min-entropy for the X and Z basis 
measurements [20, 40]. Finding higher-dimensional ana¬ 
logues of the Breidbart basis motivated Amburg et al. 
[21] to investigate MUB-balanced states. We shall see 
how qudit magic states obey many of the same desirable 
properties of the Breidbart basis. 
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III. RESULTS 

Our results are broken up into two subsections, the 
first of which concerns nonlocality while the second con¬ 
cerns entropic uncertainty. In the first subsection, we 
find a simple expression for the eigenbasis of the CHSH- 
Bell operator (13) which, after the application of a re¬ 
sult from number theory, allows us to place good bounds 
on the quantum value of the restricted (to Pauli mea¬ 
surements) CHSH game. In the second subsection, we 
examine magic states in terms of their ability to mini¬ 
mize the total (or average) entropy across a number of 
different measurement bases. For this purpose we use 
collision entropy and min-entropy, both of which have op¬ 
erational significance in cryptographic settings. Among 
the subset of qudit states that we call equatorial states, 
magic states always minimize the collision entropy, and 
for small prime dimension they also simultaneously min¬ 
imize the min-entropy. 

A. Optimal strategy for CHSH game and bounds 
on the achievable value 

Rewrite the Bell operator B* of Eq. (12) in the Weyl- 
Heisenberg basis using Eq. (14) 

B* = - Y. (20) 

and note that conjugation by a Clifford gate has the 
effect 

where 

C := CSUMi2 = Y I^X^I ® ^ ^ 

k 

is the generalized CNOT (Controlled-NOT) operator. 
Under this unitary operation the Bell operator becomes 

C^B*C=- Y (21) 

^ TiGZ*,x,sGZp 

= p5*®|0)(0| (22) 

with 

5* = - ^ (a;*(*+5)U(i|,))" (23) 

^ n^Z*,s^Zp 

and 5 = ^ ^ (a;*(*+5)U(i|,))- (24) 

so that S = S* + plp and Aniax)^) = Ai„ax(‘5*) -I-p as 
before. By the definition of stabilizer projectors in Eq. (2) 
we have 

5 = (25) 

B 


i.e. iS is a sum of projectors, one from each of the 
non-computational bases. The final step is to note that 
Amax(^*) = Amax(p5* (g) |0)(0|) = Amax(p5*) SO that 

Amax(^) = pAinax(5). (26) 

Overall, we find that the maximizing eigenvector of the 
Bell operator B gives an expectation that is the same as 
maximizing a simple single-qudit operator S. Moreover, 
iS is a sum of p stabilizer projectors, one from each of the 
non-computational stabilizer bases. We will show that 
S is maximized by a magic state |/o,6,c) as in Eq. (4). 
It follows that the optimum quantum strategy for our 
restricted CHSH game is for Alice and Bob to measure a 
shared state \Ja,b,c) of Eq. (7) since 

{Ja,b,c\B*\Ja,b.c) ={fa,b,c, 0\C^ B* C\fa,b,c, 0) 

=p{fa,b,c,0\iS*(^mO\)\fa,b,c,0) 

=p{faAc\S*\faAc) 

=P^ma.(S*) 


p 

P(l + "^) 

4p 

/g\QM 

\*^/ max 

,^xLHV 

3 

6.4641 

- 

6.4115 

6 

5 

13.9443 

20. 

13.0902 

12 

7 

22.8745 

28. 

19.4112 

19 

11 

44.1662 

44. 

34.6464 

37 

13 

56.2666 

52. 

48.3481 

47 

17 

82.9697 

68. 

55.1022 

> 66 

19 

97.4602 

76. 

72.6084 

> 79 

23 

128.508 

92. 

74.8954 

> 99 

29 

179.785 

116. 

104.819 

> 135 


TABLE I. Comparison of the bounds imposed by Information 
Causality (first column) and Weil (second column) with the 
optimal quantum and local hidden variable strategies (third 
and fourth columns). The LHV lower bounds are reproduced 
from [10] whereas the operator norm (13) is attained by 
measuring the maximally entangled states \Ja,b,c)- Because 
we have restricted ourselves to Pauli measurements, the quan¬ 
tity (13) is less than what is achievable in the unrestricted 
case. 

In Sec. IV we prove that the magic state \fa,b,c), with 
a = —1/12 and 6 = —1/8 is a maximizing eigenvector for 
the single-qudit operator 5 for p > 3. Moreover, we show 
that all magic states are balanced with respect to Alice 
and Bob’s measurements projectors given in Eq. (2) i.e., 

{KV’sl/a.b.c)!, € Zp} independent of basis B. (27) 

Our proof provides an explicit expression showing ex¬ 
actly how probabilities \{'>jj^\fa,b,c)\^ are permuted be¬ 
tween different bases. For the magic state |/-i/i 2 ,-i/ 8 ,c) 
under consideration, the overlap with vector Vb in basis 
B is constant whenever 

rB = -s(B + l). 


(28) 
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Combining this with Eq. (25) we find 

Amax(t^) = Hi’s ^ ^^^l/-l/12,-l/8,c)P (29) 

B 

= p|(+l/-l/12,-l/8,c)p (30) 

All overlaps between (non-computational) MUB vectors 
and magic states take the same form viz. 

|(^^|/aAc)|" = + + ( 31 ) 

^ k 

The values of these overlaps are constrained by the Weil 
bound [41, 42], which says that a polynomial / of degree 
n over a prime field Zp satisfies 

|^a2/(-)|<(n-l)Vp (32) 


provided p does not divide n. Applied to Eq. (30) the 
above bound implies 

(e)QM = A^a,,(5) < 4p. (33) 

Eor p > 7 this bound is more restrictive than the Informa¬ 
tion Causality bound. In Table I we compare this bound 

with XmaxiS) = (S)QM and with the bound p{l+'^^ 
imposed by Information Causality. 

The qutrit case must be handled individually - terms 
like —1/12 are ill-defined and the Weil bound is not ap¬ 
plicable. One can explicitly check that the optimal quan¬ 
tum state for Alice and Bob to measure is what we have 
come to expect i.e. a maximally entangled state that is 
isomorphic to a magic state, 

(^)S^x = (Ji.i.ol^l Jr.i.o) = 3V3cos(Vl8) = 6.4115. 

(34) 

Liang et al. [10] have shown that this is truly the quan¬ 
tum maximum; even allowing for POVMs this value can¬ 
not be surpassed. 



FIG. 2. Total min-entropy across 3 different bases: This plot 
shows contours of the total uncertainty for equatorial qutrit 
states parameterized as |^eq(/>)) = {1, with 0 < 

x,y < 9 G M. and ^ The minimum uncertainty states 

occur at the small white dots, which are integer values of x, y 
such that a; = 2a -I- 66 -I- 3c mod 9 and i/ = o -I- 66 + 6 c mod 9 
with a G Z 2 and 6 , c G Z 3 . In other words, magic states \fa,b,c) 
as defined in Eq. ( 6 ) are minimum uncertainty states amongst 
the set of equatorial states. This optimality of magic states 
\fa,b,c} holds numerically for prime dimensions p G {2, 3, 5, 7} 
but not for p = 11 or p = 13. 


The total collision-entropic uncertainty across p non- 
computational Pauli bases is bounded from below by 

Yj >-Plog I —^ I (36) 


B. Minimum uncertainty equatorial states 


Define a generalized equatorial state in terms of a vec¬ 
tor of phases (j) = (c/o := 0, (/i,..., ipp-i): 




1 

Vp 


p-1 


(35) 


This state has p — 1 free real parameters, which is exactly 
half the number of parameters that specify an arbitrary 
pure state. This is already a very important class of 
states for cryptographic protocols e.g., optimal cloning 
procedures are known for states of this form [43, 44]. 


since \ {'^eqi4‘)\V)\'^ = 1/p and we already know from 
Eq. (18) that the sum over all p -I- 1 bases is identi¬ 
cally 2 for any pure state. Because of the balanced¬ 
ness property of \fa,b,c) proven in Sec. IV we have that 
I (/a,b.civ's) 1^ independent of B and therefore the 

above inequality is saturated. In summary, of all equa¬ 
torial states (35) it turns out that magic states minimize 
the total collision entropy and saturate the above en- 
tropic uncertainty relation. 

Using magic states, we will not generally be able to 
saturate the lower bound of Eq. (19) on the total min- 
entropy across p mutually unbiased measurements. How¬ 
ever, for small prime dimensions, it is feasible to check 
numerically whether magic states minimize the total min- 
entropy when we restrict to the manifold of equatorial 
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states defined in Eq. (35). To that end define 

<?min =argmin ^ (<?))) . (37) 

$ BeZp 

For qutrits, see Fig. 2, we find that 

(fmin = {0,4’1,4>2) (38) 

27r 

= — (0,2a + 66 + 3c, a + 66 + 6c) 

whereas for primes p = 5 or p = 7 we find 

Stt 

^min = {■. ■ ■ ■) = — {ak^ + bk^ + ck) . (39) 

p 

This set of phases is exactly what defines the magic states 
\fa,b,c) in (4) and (6) In some cases the total min-entropy 
they achieve is quite close to the bound that applies to 
generic states, e.g, 

p = 3 : ^ i7^iJV'e,(<?r„i„)) = 1.468 (40) 

BGZp 

(where we used log 2 ) which is not much greater than the 
lower bound 


-3 log 2 



1.4324 


(41) 


Similarly for p = 5 we get 4.667 > 4.2113 and for p = 7 
(see Table II) we get 7.871 > 7.693. For p = 11 we 
find a counterexample to the optimality of \fa,b,c) i.e., we 
find J2Bei.p ^min\fa,b,c} = 19.8465 > 15.994 but there is 
another non-magic equatorial state that achieves 17.7606. 
For p = 13 a magic state of the form \fi^b,c) achieves 
23.471 > 20.6267 but this can also be beaten by a non 
magic equatorial state that achieves a total min-entropy 
of 23.1336. 

Minimum uncertainty states are defined relative to a 
particular notion of entropy and we have chosen two in¬ 
stances that appear to have both operational and geo¬ 
metrical interpretations. Besides the collision entropy H 2 
and the min-entropy Ha=ao, the general Renyi entropy 
Ha is well-defined for a one parameter family a > 0 , so it 
could be the case that magic states also minimize uncer¬ 
tainty for other values of a. Minimum uncertainty states 
can also be defined relative to the Shannon entropy (see 
e.g. [45]) but the results so obtained are qualitatively 
very different. 


IV. THEOREMS & PROOES 

This section contains the two key mathematical results 
and their proofs. Many of their implications have already 
been used in preceding sections. The definitions and no¬ 
tation were established in Sec. IIA 


Theorem IV.1 The magic state \fa,b,c), with a = 
— 1/12 and 6 = —1/8 is a maximizing eigenvector for 
the single-qudit operator S (25) for all prime dimensions 
p > 3. 

We claim that the eigenbasis that diagonalizes 
S is given by the set of orthogonal magic states 
{|/_i/i 2 ,-i/ 8 ,c), c G Zp}. The proof is obtained by con¬ 
structing a unitary U with these states as columns and 
showing that conjugating S with W produces a diagonal 
matrix D. 


Claim: D = U^SU = 


Ai(5) 0 0 

0 0 
0 0 Xp{S) 


where Ai(5) are the eigenvalues in no particular order. 
Perform matrix multiplication to find the explicit matrix 
elements dj^k of our putative diagonal matrix, D, i.e., 

D = J2dj,k\j){k\ = U^SU = 


j,k 


( T 2 + ^ 


^a;i2 + ®s +®'*|p)(/i| 1 >5 I ' 12 ^ 8 

g,h j \l,m 

Next, insert 


which is a consequence of Eq. (25) and Eq. (3). After 
some tidying and relabelling we find 

^j,k 

A,B,Ce’Zp 

Perform the linear substitution 


{A, B,C)^ {A + C,B - Ajl - 1/4, i - C) 

A,B,Cei,p 

which satisfies 

d'jAAj ~ 0 

because 

= pSk,j 

A&p 

Hence D is diagonal, as claimed. ■ 

Theorem IV.2 All magic states \fa,b,c) o-re MUB- 
balanced with respect to the mutually unbiased measure¬ 
ments {A, AZ,..., Moreover, the Vs-th co¬ 

efficient in the B-th basis is given by the Vo-th coef¬ 
ficient in the X measurement basis when Vb = Vq -\- 
lk{B^-^Bb) 




First define the quantity Ta^b,c as 


8 

while the maximizing eigenstate takes the form 


TaAc=-\Y. 


ak^+bk^+ck | 


P 


and note that, after a linear substitution k i—>■ k+r (where 
r is some fixed element r G Zp) that does not change the 
value of the sum, we find 


Ta,b,c = T, 


a,b-\-3ar,c-\-2br-\-3ar^ • 


(42) 


Take the probability cvo,o associated with an arbitrary 
vector Vo in the first measurement basis, B = 0: 


CVo,0 =\{4’o°\fa,b,c)\ 


p . 


afc^+&fc^ + (c+Vo)fc|2 


'^a,b,c-\-Vo 

Now consider a vector Vb in any other basis B 

CVb,B =\{'>P^^\fa,b,c)? 


(43) 


=T 

— n. 


^ a ,b—B { 2 ,c-\-Vb 


(44) 


Then, with the aid of Eq. (42), we find the condition that 
makes Eq. (43) and Eq. (44) equal: 



( 1/3 

0.7124 

0.7124 

0.0859 \ 

/l,l,o)(/l,l,0 

1/3 

0.2017 

0.2017 

0.7124 


E/3 

0.0859 

0.0859 

0.2017 / 


V. SYMMETRIES, CYCLERS AND SATO-TATE 

In this section we outline some symmetries of the states 
and operators that we have previously used. We then 
mention a connection between our work and that of Am- 
burg et al. [21] . 

To begin with, there is an obvious asymmetry associ¬ 
ated with magic states \fa,b,c)', they treat the computa¬ 
tional basis differently to the remaining Weyl-Heisenberg 
(Pauli) bases. There is nothing special about the compu¬ 
tational Z basis; in fact we can exhaustively create all all 
p‘^{p^ — 1) single-qudit |/)-type magic states by applying 
the Clifford unitary that maps Z i-G X followed by (mul¬ 
tiple applications) of the Clifford that maps X XZ. 
Any of these magic states will be balanced with respect 
to all but one Pauli measurement bases. Conversely, we 
lose nothing by restricting our analysis to the equatorial 
class of magic states that we have so far considered. 

The b and c components of magic states are modified 
by the action of Pauli operators. 


cVo,o = cvb.b, '^BGZp (45) 

^Vb = Vo + ^{B^- 4S5) ■ (46) 

A geometrical argument for the balancedness of magic 
states was implicit in recent work of Blanchfield [46] (see 
the proof of Theorem 1 therein) although that study 
concerned relationships between different MUB construc¬ 
tions. Not only have we proven that magic states \fa,b,c) 
are MUB-balanced, we have shown exactly which per¬ 
mutation of probabilities occurs when moving between 
bases. 

This theorem applied to the magic state |/-i/i 2 ,-i/ 8 ,c) 
that maximizes S gives 




which makes intuitive sense given that 

= (48) 

B 

Using the MUB decomposition introduced in the 
Sec. lie, the single-qutrit operator S looks like 


^ S GG 


I 5/3 
I 2/3 
I 2/3 


5/3 

2/3 

2/3 


2/3 

5/3 

2/3 


|/a,h,c) — ^{x\z)\fa,0,o)j 




(49) 


whereas changing the value of a in \fa,b,c) requires at 
least a Clifford operation. In fact, in [29] it was shown 
that magic states \fa,b,c) and \fa',b',c') with different val¬ 
ues a a' are connected via a Clifford if and only if these 
values lie in the same equivalence class of cubic residues 
modulo p. The noticeable effect of this is that in prime 
dimensions p = 1 mod 3 there are 3 equivalence classes of 
magic state that have different operational and geomet¬ 
rical features. In Table II we examine seven-dimensional 
magic states in terms of different operationally relevant 
quantities and see interesting differences arising from the 
inequivalent classes. 

One of the main approaches to finding minimum un¬ 
certainty states has been to find so-called MUB cyclers 
- a single unitary that, when repeatedly applied to a ba¬ 
sis state, maps basis vectors from one basis to the next, 
eventually wrapping around cyclically to return to the 
original basis. Eigenstates of these unitaries are often 
minimum-uncertainty states. These unitaries are known 
to not exist in prime dimensions, so the unitary we define 
below is the best possible alternative - it cycles through 
all but one of the bases. Define a Clifford gate Ca,b,c as 
follows 


Ca 


= M„,_b^rXMl 


b.c") 


(50) 
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then it follows from the definitions that Ca,b,c\fa,b,c) = 
\fa,b,c)- Consider C := C_i/i2 _i/8,c, then 

cn<) = iV'r’-5), K = i(-i + i) (51) 


so that for p = 5 we find that repeated application of C 
takes us through the basis indices (0, 2,4,1, 3,0, 2,4,...). 
It is fairly straightforward to prove the equivalent claim 


CIV-E) = IV'EE) VseZp, 

2 2 

(52) 

using the following argument. 


MAMtjV^r j) = 

(53) 


-l)(c+U)|fc) 


(54) 

= l/o 

(55) 

= lV'!4f) = lV'!4). 

(56) 


2 2 


Hence, magic states are eigenstates of MUB-cycling Clif¬ 
ford unitaries that cycle through p out of p -I- 1 bases. 

For prime dimensions p > 3 magic states (as defined 
in Eq. (4)) are equal weighted superpositions of roots of 
unity with a cubic polynomial in the exponent. Exponen¬ 
tial sums are commonplace in number theory and many 
results are known. The distribution 


&a,c = 




,ak'^-\-ck 


2Vp 


a e z* c e Zp 


(57) 


is real and lies in the range [—1,1] [47, 48]. The behaviour 
of Eq. (57) for different values of a, c and p is covered by 
the Sato-Tate conjecture (see [47] Sec. 1.4). Hence, magic 
states l/a,o.c) obey the same semicircular distribution of 
values -^ < (+l/a,o,c) < observed by Amburg et 
al. [21, 22, 49] for their states j/3) that are balanced 
with respect to all p -I- 1 mutually unbiased bases. Also, 
6a,c becomes increasingly equidistributed as p grows to 
infinity [48] meaning that the Weil bound Eq. (33) is 
saturated in this limit. 


VI. DISCUSSION AND CONCLUSIONS 

Characterizing highly non-stabilizer states and their 
operational capabilities is important. Stabilizer states 
appear naturally within the context of fault-tolerant 
quantum computation and also as basis vectors for com¬ 
plete sets of mutually unbiased bases. Here we have 
picked out a natural family of highly non-stabilizer states 
and showed their optimality in terms of both nonlocality 
and minimizing entropic uncertainty, where we have re¬ 
stricted to Pauli measurements in both scenarios. Find¬ 
ing concise optimal solutions to problems like these is of 
independent interest. 


a 

Wmin 

Mana 

/ j fi ^ min 

1 

-0.027692 

0.814835 

7.87055 

2 

-0.089915 

0.814835 

12.3287 

3 

-0.034531 

0.896212 

9.35125 

4 

-0.034531 

0.896212 

9.35125 

5 

-0.089915 

0.814835 

12.3287 

6 

-0.027692 

0.814835 

7.87055 


TABLE 11. Different characterizations of seven-dimensional 
magic states : Magic states have surprising additional struc¬ 
ture in prime dimensions of the form p = 1 mod 3 [29], 
which opens the door for operational inequivalence of magic 
states \fa,b,c) depending on the value of the parameter a £ 
{l,2,...,p— 1}. The nonzero elements of Zy can be par¬ 
titioned into equivalence classes of cubic residues and their 
cosets i.e., Zy = {1,6} U {2,5} U {3,4}. Depending on 
which of the operationally relevant quantities {Wmin, Mana 
or 50 b Tlmin) we focus on, we find a different choice of opti¬ 
mal equivalence class. The quantity Wmin denotes the min¬ 
imal value of the Wigner function of \fa,b,c} in phase space 
[17], while Mana effectively describes the magnitude of the 
sum of all the negative entries [18]. The last column de¬ 
scribes the total min-entropy of \fa,b,c.) across all measure¬ 
ment bases {X, XZ,..., X Note that lower bound 
for arbitrary (not necessarily equatorial) states is 7.693 = 
—71og[(l -I- (7 — l)/%/7)/7] so the magic state performs well. 


Whenever Alice and Bob are restricted to using qudit 
Pauli measurements we cannot ever observe nonlocality 
by measuring a stabilizer state. This is due to the exis¬ 
tence of a local, noncontextual hidden variable model for 
non-negatively represented states in a particular discrete 
Wigner function [16, 17, 51]. By showing that magic 
states give the maximal violation relative to a stabilizer 
CHSH scenario, this gives an operational interpretation 
to the concept of being highly non-stabilizer. Similarly, 
a stabilizer state is a maximally uncertain state with re¬ 
spect to a set of mutually unbiased Pauli bases. By show¬ 
ing that magic states minimize the entropic uncertainty, 
this gives another operational characterization of their 
highly non-stabilizer character. 

One could argue the case for a number of different 
ways of quantifying non-stabilizerness. The correct met¬ 
ric will most likely depend on the exact nature of the 
task for which this state or operation is a resource. In 
[18] two measures were put forward that were tailored 
toward fault-tolerant quantum computation via magic 
state distillation. The relative entropy between the test 
state and the closest positively represented state, as 
well as the sum-negativity (the sum of all the nega¬ 
tive quasi-probabilities in the discrete Wigner function) 
of the state were both highlighted. It is worth noting 
the disagreement in Table H between the state l/ 3 ,fc,c) 
that Mana picks out as opposed to the state l/i,6,c) that 
min-entropy picks out. Other studies have characterized 
non-stabilizer states and operations in terms of convex 
geometry [55] and in terms of frame potentials [56]. 
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For qubits, the Tsirelson bound for the CHSH sce¬ 
nario can be saturated using Pauli measurements applied 
to a particular maximally entangled state. The maxi¬ 
mally nonlocal state |Ti,o,o) (H) can be understood as 
the image of a “pi-over-eight” gate applied to one half 
of a Bell pair |<f>). It is intuitively pleasing that this 
connection between optimality and the Clifford hierar¬ 
chy continues for all higher prime values of p, as we have 
shown here. Less pleasing is that the quantum advan¬ 
tage diminishes with increasing dimension, and seems 
to disappear entirely for p > 11 (we highlight as an in¬ 
teresting open question whether any alternate choice of 
Pauli measurements (14) for Alice and Bob could result in 
Amax(^) > (^)max = 37). For p < 11 it should be borne 
in mind that [13] showed a clear operational advantage 
for restricting to stabilizer measurements when statisti¬ 
cal significance, rather than the amount of Bell inequality 
violation, is adopted as the figure of merit. One possible 
practical application of our results is the benchmarking 
of non-stabilizer resources for fault-tolerant computation. 
Any imperfect \Ja,b,c) (or equivalently Ma^b,c) that still 
allows for violation of the qudit Bell-CHSH inequality is 
also suitable for promoting fault-tolerant Clifford gates to 
universal quantum computation via magic state distilla¬ 
tion [50] . This is a sufficient condition but not necessary, 
as the distillation routines of [26] can also enable univer¬ 
sality with imperfect \Ja,b,c} that have lost the ability to 
violate the Bell inequalities [9] that we used here. 

The inter-relationship between entanglement, nonlo¬ 
cality, steering, complementarity and uncertainty is a fas¬ 
cinating and ongoing research program. Oppenheim and 
Wehner [53] have shown that the presence of uncertainty 
limits the amount of nonlocality, whereas Tomamichel 
and Hanggi [54] have shown that uncertainty is necessary 
in order to observe nonlocality. In the qubit case, the 
similarity between the Information Causality bound of 


Eq. (17) and the entropic uncertainty bound of Eq. (19) 
arises because of these links between nonlocality and un¬ 
certainty. The fact that the qudit expressions also have 
the same form suggests that this connection may hold 
generally, although this is complicated by the fact that 
the qudit CHSH game is no longer an XOR game [11]. 

There are some obvious avenues for further investiga¬ 
tion. Magic states \fa,b,c) reflect the strncture of the 
third level of the Clifford hierarchy, and take a simple 
exponential sum form for higher dimensions. Both of 
these facets are amenable to further analysis using group- 
theoretical and number-theoretical ideas. Moreover, they 
both suggest natural generalizations for classes of states 
that may prove useful, analagous to what we have proven 
here. The Clifford hierarchy is well-defined for multi¬ 
ple particles (e.g. controlled-Clifford gates are elements 
of the third level) but for single-particle diagonal gates 
the relationship between hierarchy level and the order 
of polynomials in the exponential sum is more transpar¬ 
ent. The Weil bound (33) suggests that states defined as 
exponential sums of higher order polynomials may some¬ 
times perform better. Finally, as with all calculations 
over prime fields, it is worth investigating how many of 
our results carry over to the case where we consider sys¬ 
tems of prime-power dimensions q = p'" hy using the field 
Fq rather than Fp. 
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